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Exact relationship between the entanglement entropies of XY and 
quantum Ising chains 
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Abstract. - We consider two prototypical quantum models, the spin- 1/2 XY chain and the 
quantum Ising chain and study their entanglement entropy, S(l, L), of blocks of £ spins in homo- 
geneous or inhomogeneous systems of length L. By using two different approaches, free-fermion 
techniques and perturbational expansion, an exact relationship between the entropies is revealed. 
Using this relation we translate known results between the two models and obtain, among others, 
the additive constant of the entropy of the critical homogeneous quantum Ising chain and the 
effective central charge of the random XY chain. 



Introduction. — Recently, we have witnessed a grow- 
ing interest in entanglement effects in quantum many- 
body systems [1]. If an isolated quantum system is di- 
vided into two parts, A and B, all information about A is 
contained in the reduced density matrix: pa = Trg |0}(0|, 
where |0) denotes the ground state of the system. To quan- 
tify the quantum entanglement between A and £>, different 
measures have been introduced; a frequently used quantity 
is the von Neumann entropy defined by: 



Sa = -Tr(p A log 2 pa) = -Tr(p e log 2 ps) 



(1) 



Contrary to the thermal entropy, the entanglement en- 
tropy is not an extensive quantity, but for a non-critical 
system it scales with the area of the interface [2] separating 
• the block (i.e. A) and the environment (£>). In the case of 
a one-dimensional (Id) infinite system and a finite block 
of length, I, the interface consists of a few points and the 
entanglement entropy generally approaches a finite value 
in the limit I — > oo. At the critical point, however, Sa 
diverges logarithmically as 



Cl, 



(2) 



where the prefactor, c, is universal in homogeneous sys- 
tems and given by the central charge of the associated 
conformal field theory [3], whereas the constant, c\, is non- 
universal. The relation in Eq.(2) has been generalized by 
Calabrese and Cardy [3] for finite systems, for different 



boundary conditions, for non-critical systems in the vicin- 
ity of the transition point and for finite temperatures. In 
higher dimensions only a few results are available, mainly 
for non-interacting fermions [4] and bosons [5,6] but also 
random quantum Ising models have been investigated [7]. 

Some of the conformal results have been tested on inte- 
grable quantum spin chains in particular on the antiferro- 
magnetic AF-chain, defined by the Hamiltonian [8]: 



H 



XY 



E 

i=l 



(JfSfS: 
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(3) 



J L+1 — 



Here the S^' v, s are spin- 1/2 operators at site i 
Si' v and the couplings Jf and Jf may be different and 
site-dependent in general. In the followings, we restrict 
ourselves to even system sizes L. If the interaction is 
isotropic on average in the sense that [In J x ] av = [In J v ] av 
holds, where [. . .] av stands for the average over the distri- 
bution of couplings, the model is critical, which manifests 
itself in the vanishing of the gap and the algebraic decay of 
correlation functions in the thermodynamic limit L — > oo. 
For the special case of the homogeneous AA-chain, i.e. for 
Jf — Jf — J both the prefactor (c = 1) and the constant 
ci in Eq.(2) have been calculated exactly [10]. 

Another basic one-dimensional quantum model is the 
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quantum Ising chain (QIC) defined by the Hamiltonian 



1 L 1 L 



(4) 



in terms of the Pauli-operators oy z = 2S*f' z at site i 
(tr^'fj = cr^' z ) and Aj and /ij are nearest neighbor cou- 
plings and the transverse fields, respectively. This model 
exhibits a quantum phase transition if [lnA] a „ = [ln/i] a „, 
i.e. for the homogeneous system with Ai = A and hi = h, 
the critical point is located at J = h [9]. The satura- 
tion value of the von Neumann entropy in the homoge- 
neous model has been found to exhibit a singularity in the 
vicinity of the critical point as Soo = | log 2 £, where £ 
denotes the correlation length diverging at criticality [3]. 
The prefactor is obtained to be c = 1/2 in an analytical 
calculation using a mapping between the reduced density 
matrix of the model and the corner transfer matrix of the 
2d classical Ising model [3, 11]. 

The entanglement entropies of the XX and the quan- 
tum Ising model have also been studied in the presence 
of quenched disorder, when the parameters of the models 
(the isotropic couplings Jf = Jf = Ji, or the Aj bonds and 
the hi transverse fields, respectively) are independent and 
identically distributed random numbers. The average en- 
tropy has been calculated analytically [12] by a strong dis- 
order renormalization group method [13], which is found 
to follow the logarithmic law in Eq.(2). In this case the 
prefactor, which is called the effective central charge and 
denoted by c e ff is obtained as c c tf(XX) = In 2 for the ran- 
dom AA-chain and c e s(I) = In 2/2 for the random critical 
QIC, respectively. 

Numerical studies of the entanglement entropy were 
performed on different XX-chains, in which the effect of a 
free boundary [14], a defect coupling [15], random [16] or 
aperiodic interactions [17] and the presence of an energy 
current in the system [18], etc. were investigated. For the 
random QIC, the location of the maxima of the entropy is 
used to define sample-dependent critical points [19]. The 
evolution of the entropy after a quench in both models is 
also the subject of recent investigations [20,21]. 

It is known for some time that the two Hamiltonians in 
Eqs.(3) and (4) can be mapped into each other through 
a canonical transformation [22-24], which is described in 
the Appendix. As a consequence, the spectrum of the 
two Hamiltonians, as well as some correlation functions of 
the two models are related, as well. One might ask the 
question, whether a similar correspondence can be found 
concerning the entanglement entropies of the two models. 
At first thought, the existence of such a relation is not ob- 
vious since the transformation of the operators in Eqs.(34) 
and (35) is nonlocal, c.f. operators in A for one model are 
expressed with operators located both in A and B. 

In this paper, we study in detail the relation between 
the entanglement entropies of the two models by two ap- 
proaches. In the first approach we calculate the entan- 
glement entropies in the free-fcrmion representation [8]. 



By this method, the Hamiltonians are first expressed in 
terms of free fermions, which requires the solution of an 
eigenvalue problem of dimensions, Lx L, then after a sec- 
ond transformation, the systems assume the form of non- 
correlated fcrmionic modes, which is obtained by solving 
an I x I eigenvalue problem. These eigenvalue problems 
are then compared for the two models. In the second ap- 
proach, the entanglement entropy is calculated by a per- 
turbation expansion, in terms of different powers of the 
coupling term connecting the two parts of the system and 
the expressions obtained for the two models are then com- 
pared. The relation between the entropies is then used to 
transfer existing results between the two models. 

Free-fermion calculation. — The main steps of the 
calculations to be carried out in this section are sum- 
marized as follows. The Hamiltonians are expressed in 
terms of fermion operators and are diagonalized. Then, 
the restricted correlation matrix is constructed and it is 
transformed to a form corresponding to non-correlated 
fcrmionic modes by a canonical transformation. Finally, 
the entanglement entropy is calculated from the eigenval- 
ues of this matrix. 

Diagonalization of the Hamiltonians. Both models 
can be expressed in terms of fermion creation and an- 
nihilation operators c\ and c,, respectively [8], which 
are obtained through the Jordan- Wigner transformation: 



exp 



Ei-1 + - 
7 a i a i 



and c+ = af 



exp 



Ei— l + - 
3 "J"! 



Ci = 



The Hamiltonian of the XY- 



chain in Eq.(3) is expressed as: 



L-l 



n X Y = J2\{(Jf-J?)ctcf +1 + (J? + JfK+C i+1 +h.C.} 
l—l 

\w{{Jt Jl)cl4 + (Jl + Jl)c+ Cl + h.c.} (5) 

where w = cxp(i7rA c ), with N c = 1 cfcj. 

Similarly one obtains for the Hamiltonian of the QIC [9] 
inEq.(4): 



Hi = X>k + 



i=l 
1 



1\ I" 



i=l 



+ ^w\ L (c+ -c L )(c+ +ci). 



(6) 



For both Hamiltonians, the N c is even in the ground state, 
thus we have w — 1. Moreover, both Hamiltonians are 
quadratic in c+ and Cj, therefore they can be diagonalized 
by standard techniques [8]. In our approach, the basic 
quantity is a 2L x 2L matrix, denoted by T, which can be 
interpreted as the transfer matrix of directed walks and 
which has been introduced in [25] for the QIC and in [24] 
for the XY-chain [26] . From the solution of the eigenvalue 
problem of T one obtains both the energies of the free- 
fermionic modes, Afe > 0, k = 1,2, ... L, and two sets of 
vectors: and \&/c, both having L components [27]. For 
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the XY-chain, T is given by [24]: 

see eq. (7) above 

and the eigenvectors contain the 
(* fc (l), tffc(l), * fc (2), * fc (2), $ fe (3), * fc (3) . . . 
* fc (L-l),* fc (L),$fc(L)). 

For the QIC the Tj matrix reads as [25]: 



components: 



T, = 



hi 

Ai 
Ai 



-w\ L \ 



\L-1 











(9) 



and the eigenvectors have the components: 
(-$ fc (l), * fc (l), -* fc (2), * fc (2), . . . , -$ fc (L), * fc (L)). 

As shown in [24] the L eigenvectors of Txy with posi- 
tive eigenvalues can be divided into two classes: 

i) For the first class of vectors, which we mark with 
odd superscripts we have $2fe-i(2i) = \^2fc-i(2i — 1) = 
0, whereas the non-zero components of the vectors are 
obtained from the eigenvalue problem of the matrix: 



( o 

Jl 



\-wJi 



J y 
J i 





JX 

J 2 



T y 



Jl-2 



-wJf\ 



J y 

J L-1 



o / 

(10) 

and the eigenvectors have the components: 
(* fc (l), * fc (2), $ fc (3), * fc (4), . . . , § k {L - 1), tt fc (L)). 

This is just the T matrix of a QIC with the Hamiltonian 
Hi (a) given in Eq.(36). Denoting the components of the 
vectors corresponding to this QIC by <f>£\i) and ^^\i) 
we have the correspondences: 



$ 2fc _ 1 (2i-l) = -$^ ) (i), * 2 fc-i(2i) 



(11) 



ii) For the eigenvectors of the second class, which are 
marked with even superscripts, the vanishing components 



are $2/0(2* — 1) = ^>2fc(2i) = 0, whereas the non-zero com- 
ponents of the vectors are obtained from the eigenvalue 
problem of the matrix: 



T (r) / 



/ 

Jf 



\-wJ 



Jf 


J 2 



T y 

J 2 



T y 

J L-2 





jx 

J L-1 



-wJ y \ 



TX 

J L-1 



J 

(12) 

and the eigenvectors have the components: 
(* fc (l), $ fc (2), * fc (3), $ fc (4), . . . , V k (L - 1), $ fc (L)). 

This is again the T matrix of a QIC with the Hamil- 
tonian Hj{t) given in Eq.(36). Denoting the components 
of the vectors corresponding to this QIC by <&^\i) and 
^>^\i) we have the correspondences: 



$ 2 fc(2i) = * 2fe (2z-l) = -$i r) (i 



(13) 



Thus we conclude that the explicit solution of Txy re- 
quires the diagonalization of two T-s of QIC-s, which is 
just equivalent to the mapping described in the Appendix. 

Correlation matrix and entanglement entropy. Next 
we consider a block of length £, consisting of spins i = 
1,2 ,...,£ and the reduced density matrix is given by: 
Pi = TrL-e |0)(0|. For free fermionic systems pe can be re- 
constructed from the restricted correlation matrix [28,29], 
G, the matrix-elements of which are given by: 



G v 



= (0|(c+-c„)( C ++c m )|0) 

L 

= -^* fe (m)*fc(n), m,n=l,2,...,£ 



(14) 



k=l 



For the XY chain using the properties of the vectors, <frfc 
and *fc in Eqs.(ll) and (13) we obtain for the matrix- 
elements: 



G2i,2j — 0, 



G2i-l,2j-l — 
■»(ct) r< _ _/^( r ) 



(15) 



where G 



denotes the matrix-element of the correlation 



matrix of the QIC with Hamiltonian Hi (a, r) . The cor- 
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relation matrix for even £ is bipartite, being composed of 
2x2 matrices 





-a 



-G 




(r) 

3,i 



(16) 



i,j = 1,2,... e/2. 

In order to obtain the von Neumann entropy, S(£,L), 
one diagonalizes pi, which is given through a canonical 
transformation: 



i 



- (v q (i) + u q (i)) c l + - (v q (i) - u q {i)) c+ 



(17) 

where the v q (i) and u q (i) are real and normalized: 

Si v q(i) = Y^i u 'q(i') = 1- ln tnc transformed basis we 
have 



(0|/x ?Mp |0) = 0, 



I Mo Mp|0> = V 



1 



(18) 



for p, g = 1,2,... I. Thus the fermionic modes are mi- 
correlated and the eigenvalues of are the products of 
(1 ± v q )/2, q = 1,2,... I. The entropy of the system is 
given by the sum of binary entropies: 



S{t,L) = -Y J 

9=1 



1 



Vq , 1 

— !og 2 



+ 



1 



2 oz 2 2 
The v q -s are the solution of the equations: 



Vq 1 

— !og 2 - 



Gu q = U q Vg, G T V q = V q U q , 



or, equivalcntly, one has 



GG T v„ 



v\v q , G T Gu, 



(19) 
(20) 

(21) 



For the XY chain with even £, the matrix GG T is com- 
posed of 2 x 2 diagonal matrices 



[( G (-)) T G^)] 







[GW(GW) 1 



(22) 



thus GG T can be written in a form which consists of two 
diagonal blocks and the eigenvalues are obtained by solv- 
ing two separate eigenvalue problems of (G^) T G^ and 
(G^) T G (T K Now, it follows from Eq.(19) that the en- 
tanglement entropy of the XY-chain in Eq.(3) is the sum 
of the entanglement entropies of the two QIC-s defined in 
Eq.(36): 

S {XY) {£, L) = {£/2,L/2) + 5 (T) (£/2, L/2) . (23) 

This relation constitutes the main result of our paper. 

One can make sure easily that the above result holds 
not only for blocks of contiguous spins but for any block 
composed of pairs of adjacent sites (2i — l,2i). In any 
other case, the matrix GG T is still block-diagonal and 
the von Neumann entropy can be written as a sum of two 
terms each of which depend exclusively on the parameters 



of one of the decoupled QIC's, however, these terms are 
no longer to be interpreted as entropies of some blocks in 
the QIC's. 

Next, we turn to the case of odd £, when the correlation 
matrix is obtained by leaving out the last row and column 
of the matrix with £ + 1, thus the structure in Eq.(16) 
is lost. The matrix G consists of non-square submatri- 
ces of size / x (I + 1) and (I + 1) x I of G< T ) and G^\ 
respectively. Consequently, the eigenvalues of GG T can 
not be expressed by those obtained from the two QIC-s. 
Then the relation in Eq.(23) is only asymptotically valid, 
as £ ^ 1 and the corrections are of the order of l/£. 

Perturbative calculation. — In this section we split 
the Hamiltonian of the system as: 



n = h a + n B + v , 



(24) 



where Ha and Hb are the Hamiltonians of the free sub- 
systems, A and B, respectively, and V is the interaction 
term, which reads for the XY-chain as: 



V XY = H(£)+H(L) 
H{€) — J g S\ Sgj^i + Jg Sg Sg, 
H(L) = JlSfSf + J y L SlSf . 



(25) 



Let us denote the eigenstates of Ha by \(pf) with energies 
Ef, and similarly for Hb the eigenstates are \tp%) with 
eigenvalues Ef. The ground state of the total system with 
Hamiltonian H can be expressed in terms of \<p^) and \(fi®) 

as l u ) = J2t J2k c (*' k)(\<pf} ® l^f))> so that, the reduced 
density matrix, pg, has the matrix-elements: 



{<rt\pt\<pf)=pt{i,j) = E«^l®(^l)0)<0(|^>(8)|^» 

k 
B 

= Y,c{i,k)c*(j,k) . (26) 



Here c(0, 0) = 1, otherwise the expansion coefficients 
c(i, k) are calculated perturbatively. In leading order we 
have 

and the higher order terms are sums of products contain- 
ing factors of the form: 

/(i , i1wl = MW»iM, (28) 

Ef + E*-E$-E* K ' 

such that Ef + El - E$ - E% ± 0. For the XY-chain 
with the interaction term in Eq.(25) we have f(i, i'\k, k') = 
fe(i, i'\k, + i'\k, k'), in which the first term is given 
by: 

Mi,i%k') = -^[Jf{<fif\Sf\<tf){tf\Sf +1 \<p%,) 



+ jy{vt\sj\rf){v B k \s v l+ M,)] (29) 
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with AE = Ef + E% - E$ - E®, and there is a similar 
expression for ,/l(i, i'\k, k'), as well. Note that the matrix- 
elements in these expressions are separated as the product 
of two matrix-elements of the surface operators in the two 
subsystems. 

In the next step we perform the same perturbation ex- 
pansion in terms of two independent QIC-s with parame- 
ters given in Eq.(36), in which case for even £ the pertur- 
bation is located at £/2 and L/2 and given by: 

H{£/2) = -^J£Vt/2< J e/2+i + ^Jf T e/2 T e/2+i 
H(L/2) = \j x L al /2 al + \j v L Tl /2 Tt . (30) 

Using the mapping in Eq.(35) we obtain H{£) = H{£/2) 
and 

L/2 L/2 

H(L) = -Jlol^l + - J v L rl /2 rf ]J a] . (31) 

3=1 3 = 1 

For the QIC calculation we denote by ff/ 2 (i,i'\k,k') the 
factor, in which we use the states labeled by i,i'\k, k! 
in Eq.(28). Then AE remains the same, as well as 
the matrix-elements are invariant, so that fg(i, i'\k, k') = 
fi/ 2 (hi'\k,k'). For the other term, fL/ 2 (hi'\k,k'), the 
only difference, that the transformed perturbation con- 
tains also the products IIj=i T j an d rij=i a j i which com- 
mute with the Hamiltonians -Hj(t) and Hi (a), respec- 
tively, and have the eigenvalues p = ±1. The excited 
states, however, which enter into the expansion of c(i, k) 
have the same parity as the ground state, p=l, thus also 
f L {i,i'\k,k')=~f L/2 {i,i'\k,k'). 

We conclude that for even I the expansion coefficients 
c(i, k) are identical for the XY model as well as for two 
independent QIC-s with parameters given in Eq. (36) in 
all order of the perturbation expansion. Consequently, the 
same is true for the elements of the reduced density matrix 
and finally for the entanglement entropy. In this way, we 
have reobtained the result already calculated in Eq.(23). 

If the size of the cell, I is odd, then the interaction term 
in Eq.(25) is transformed as: 

H {£) = ^t (J (t+\)/2 + 4^ T (£+l)/2 ( 32 ) 

and similarly for H(L). These cannot be written in a 
separated form in terms of the a and r operators, thus 
Eq.(23) is no longer valid. 

Discussion. — In this paper we have derived an exact 
relation in Eq.(23) between the entanglement entropy of 
the XY-chain and that of the QIC. This relation is valid for 
a finite block of even size and holds also for inhomogeneous 
couplings. Since the derivation is based on a mapping 
between the two correlation matrices similar relations can 
be obtained for another measures of the entanglement, 
such as the Renyi entropy or the concurrence. 



Before discussing the simple consequences of the rela- 
tion in Eq.(23) we begin with the comb entanglement [30] 
of the XY-chain, when the block A consists of t < L/2 
spins, which occupy sites having the same parity. Then, 
according to Eq.(15) all elements of the matrix G are 
zero and the entropy is S(£, L) = £, so that the block 
is maximally entangled with the environment. In case of 
strongly disordered XX chains, which possess asymptoti- 
cally a random singlet ground state [13], this finding can 
be obtained directly since singlet bonds form exclusively 
between spins at sites with different parities. As can be 
seen, the singlet-state approximation happens to give an 
exact result for the comb entanglement. 

First we consider a homogeneous XY chain with J x = 
(I+7) and J y = (I+7), so that in the equivalent decoupled 
QIC-s we have A = 1 and h = (1 — 7)/(l + 7) (h = 
(l + 7)/(l — 7)) for the a (r) chain. In the thermodynamic 
limit according to Eq.(23) we have: 

S <-»( 7 ).5»»(,^i T 2) +S ''>(^i±2). 

(33) 

which is indeed satisfied with the known exact results 
[10,11]. For the XX-chain, which corresponds to 7 = 0, 
the entropy is just the double of the entropy of the criti- 
cal QIC. Consequently, the central charges are related as 
c(XX) = 2c(I), which is indeed observed in the exact 
calculations [3,10,11]. For the non-universal constant in 
Eq.(2) we have the relation a(XX) = 2[-^ + Ci(J)]. 
Since c\(XX) is also exactly known [10], in this way we 
have obtained the exact value of c\ (I) , which has not been 
known in the literature. 

Our next remark concerns the entropy profile of a fi- 
nite and open XX-chain, i.e. with Jl = 0, which has 
been recently calculated [14] and a staggered behavior is 
observed depending on the parity of £. According to our 
result for even £ the mapping with the QIC is perfect, 
whereas for odd £ there are finite-size corrections, which 
scale with l/£. This result is in complete agreement with 
the observed numerical findings. 

For disordered chains the relation in Eq.(23) is valid 
for each (set of) samples, consequently the average en- 
tropies are also related by a factor of two. If the disorder 
is in the XX- form, i.e. Jf = J? = Ji, then the effec- 
tive central charges satisfy: c c s(XX) = 2c e g(I), which is 
also in agreement with the results obtained by the strong 
disorder renormalization group method [12]. We obtain, 
however, new results, if the disorder is XY-type, i.e. gen- 
erally Jf ^ Jf. According to our mapping the effective 
central charge of random critical XY-chains is the same as 
for random XX-chains and given by: c e s(XY) — In 2. 

The effect of another type of inhomogeneitics on the 
entanglement entropy have also been considered. In this 
respect we mention numerical studies of the XX-chain hav- 
ing one or two defects, connecting the two subsystems [15], 
or quasi-periodic, or more generally aperiodic modulation 
of the couplings [17]. These perturbations result in varying 
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effective central charges in both models which are however 
related through Eq.(23). 

Note added in proof: After submitting this work a pa- 
per by Cardy et al has appeared [31] in which c\(I) is 
calculated by different methods. 
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Appendix: Mapping of the XY chain to two de- 
coupled QIC-s. — This type of mapping in the thermo- 
dynamic limit is presented in [22] and applied for random 
chains in [23]. For finite chains it is described in [24]. 

Let us define two sets of Pauli operators, cr^' z and ff ,z , 
i = 1,2,..., L/2 through the spin-1/2 operators Sj' v , with 
j = l,2,...,Lby: 



of=l[(2S?), <=4S 2 V^ 

2i-l 



(34) 



The inverse transformations are given by: 

i— 1 i 



n 



25|i — r- 



n 



(7 J 



(35) 



In terms of these Pauli operators the Hamiltonian operator 
of the XY-chain with L spins in Eq.(3) can be written as 
the sum of two decoupled quantum Ising chains with L/2 
sites: 



Wi(r) 



L/2 



L/2 



- 2 X] J 2*«+l - 2 X/ J 2i-l°f 
L/2 L/2 

J V T X T X — T X T Z 



(36) 



i=l 



Here, in the last step, we have made a gauge transforma- 
tion to change the sign of the right-hand side of the last 
two equations. Note that although [Hi(cr), Hi(t)] = 0, 
the two chains are not independent since e.g. of and rf 
do not commute with each other. 
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